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ABSTRACT 

We derive [K/Fe] abundance ratios for 119 stars in the globular cluster NGC 2808, all of them having O, Na, 

Mg and A1 abundances homogeneously measured in previous works. We detect an intrinsic star-to-star spread 
in the Potassium abundance. Moreover [K/Fe] abundance ratios display statistically significant correlations 
with [Na/Fe] and [Al/Fe], and anti-correlations with [O/Fe] and [Mg/Fe], All the four Mg deficient stars 
([Mg/Fe] <0.0) discovered so far in NGC 2808 are enriched in K by ~0.3 dex with respect to those with normal 
[Mg/Fe], NGC 2808 is the second globular cluster, after NGC 2419, where a clear Mg-K anti-correlation is 
detected, albeit of weaker amplitude. The simultaneous correlation/anti-correlation of [K/Fe] with all the light 
elements usually involved in the chemical anomalies observed in globular cluster stars, strongly support the 
idea that these abundance patterns are due to the same self-enrichment mechanism that produces Na-O and 
Mg-Al anti-correlations. This finding suggests that detectable spreads in K abundances may be typical in the 
massive globular clusters where the self-enrichment processes are observed to produce their most extreme 
manifestations. 

Subject headings: stars: abundances — techniques: spectroscopic — globular clusters: individual (NGC2808) 


1. INTRODUCTION 

In the last decade, accurate photometric and spectroscopic 
investigations have changed our concept of globular clusters 
(GCs), which were considered for a long time as aggregates of 
coeval stars, born with the same initial chemical composition. 
This traditional paradigm is still valid in terms of the iron and 
iron-peak elements abundance, while large variations in light 
elements (O, Na, Mg and Al) were revealed in all the GCs 
studied so far, in our Galaxy (see e.g. ICarretta et al.l 120091: 
Gratton, Carretta & Bragaglia 2012) and in the Local Group 
dMucciarelli et all20C)9t) 

Chemical inhomogeneities are usually explained by invok¬ 
ing the existence of (at least) a second generation of stars 
formed within the first ~ 100 Myr of the cluster life, from 
intra-cluster medium polluted by the first generation stars. 
Intermediate mass Asymptotic Giant Branch (AGB) stars 
dD’Ercole et al .1120081) . Fast Rotating Massive Stars (FRMS 
iDecressin et al.(l2007l) and binary stars (De Mink et akll2009h 
were proposed as the most likely polluters, that should signif¬ 
icantly affect the abundance of light elements (e.g. enhancing 
Na, Al and He and depleting Mg and O) while leaving the 
iron abundance unaffected. However, all these models face 
serious problems and we are still far from a full understand¬ 
ing of the processes underlying the presence of multiple pop¬ 
ulations in GCsJsee e.g. the alternative scenario proposed by 
iBastian et al.ll2013l) . 

1 Based on data obtained at the ESO Very Large Telescope under the 
programs 072.D-0507 and 091.D-0329. 
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Potassium is a new entry among the chemical anoma¬ 
lies in GCs. From the analysis of DEIMOS@Keck low- 
resolution spect ra of 49 giant stars in the remote, massive 
GC NGC 2419. lMucciarelli et alJd20ll Mul2 hereafter) dis¬ 
covered that its stars span an unusually large range in K 
abundances, from solar values up to [K/Fe]~ +2 dex. Also, 
[K/Fe] anti-correlate with [Mg/Fe] that also covers a huge 
range, from values compatible with the ct-enhancement ob¬ 
served in other ancient GCs down to [Mg/Fe] ~ -1 dex. About 
40% of the stars in NGC 2419 show sub-solar [Mg/Fe] abun- 
da nce ratios and high K abundances. This were confirmed 
bv ICohen & Kirbvl (120121) from the analysis of HIRES @Keck 
high-resolution spectra of 13 giant stars of NGC 2419. Note 
that most of the GC stars have [Mg/Fe] abundance ratios 
between +0.3 and +0.6 dex (thus, with variations generally 
smaller than those observed in the other light elements), and 
only a handful of GC stars with -0.3<[Mg/Fe]< +0.3 are 
known (I Carretta et al]l2009l) . NGC 2419 represents the only 
exception observed so far. 

Currently, a solid interpretation of the origin of such a 
Mg-K anti-correlation is still lacking. I Ventura et all (12012 L 
V12 hereafter) proposed a theoretical model where the Mg- 
poor/K-rich stellar population is an extreme population di¬ 
rectly formed from the AGB and super-AGB ejecta, support¬ 
ing the first claim by Mu 12 that the spread in K is ascrib- 
able to the same self-enrichment process able to produce the 
other chemical anomalies. However, a proper modeling of the 
observed pattern requires a fine tuning in the reaction cross- 
sections and burning temperatures. 

ICarretta et all (|201 3l) analysed small samples of stars in dif¬ 
ferent evolutionary stages (turnoff, sub-giant and giant stars) 
in seven GCs, finding no evidence of intrinsic scatter in their 
K content. All the studied stars have Mg and K abundances 
that well match the values of the normal stellar population of 
NGC 2419. 

In order to unveil the possible role played by K in the frame¬ 
work of the GC self-enrichment process, and to understand 
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whether the Mg-K anti-correlation in NGC 2419 is a singular 
event among the chemical anomalies of the GCs, we started a 
project aimed at deriving the K abundance in different GCs. 
This Letter is focused on the massive GC NGC 2808, one 
of the clusters with the most extended Na-O anti-correlation 
(see e.g. I Carretta et all1200 61) and h arboring 4 Mg-poor stars 
(ICarretta et al J2009l : ICarrettal20 1 4l) . making it an ideal cluster 
where K abundance spreads can be searched for. 

2. OBSERVATIONS 

The observations were performed with the spectro¬ 
scopic facility FLAMES dPasquini et ali |2000|) in the 
UVES+MEDUSA combined mode. We adopted the 860 Red 
Arm CD4 UVES set-up, with a spectral coverage of 6600- 
10600 A and a resolution of ~45000, and the HR18 GI¬ 
RAFFE grating, covering from 7648 to 7889 A and with a 
spectral resolution of 18400. These set-ups were chosen in 
order to sample the K I resonance line at 7699 A . 

We used two fibre allocation configurations to observe 
our targets, that are 119 Red Giant Branch (RGB) stars 
whose membership to NGC 2808 was already established, 
and whose Fe, O, Na, Mg, and Al abundances we re homo¬ 
geneously estimated by Carr etta et all {2006, 2009, C06 and 
C09 hereafter, respectively). The goal of our project is to build 
on their analysis by adding the abundance of Potassium. In 
the following all the adopted abundance ratios are from these 
works, except for [K/Fe]. 

For each fibre configuration two exposures of 1300 sec each 
were secured, in order to reach a S/N ratio per pixel around the 
K line of ~ 100 for the faintest targets (V~ 15.4) and of ~200 
for the brightest ones (V~13.8). All the spectra were reduced 
with the UVES-FLAMES and GIRAFFE ESO pipelines, in¬ 
cluding bias subtraction, flat-fielding, wavelength calibration, 
spectral extraction and order merging (only for UVES spec¬ 
tra). 

3. K ABUNDANCES 

The K abundances of the targets were derived from the 
measure of the K I line at 7699 A with the package GALA 
(iMucciarelli et al.l 120131) . The line equiv alent widths (EWs) 
were measured with the code DA OSPEC (IStetson & Panciiiol 
f2008l) . through the wrapper 4DAO (lMucciarellil2013i) . Eleven 
stars were observed both with UVES and GIRAFFE spectra: 
we found an average difference between their K line EWs of 
£Wuves - £Wgir = 2.5 ±1.2 mA (a = 4.0 mA), corresponding 
to a typical variation smaller than 2% in EW (~0.03-0.04 dex 
in K abundance). 

Effective temperatures (7/n) and surface gravities (log g) 
are from C06. The analysis of the entire FLAMES sample 
of 4 GCs secured within this project (Mucciarelli et al. in 
prep.) has revealed that the use of the micro-turbulent veloc¬ 
ities (vturb) derived by C06 leads to an unexpected trend be¬ 
tween K abundances and v tur b, probably due to the small num¬ 
ber of Fe I lines used to infer v tul -b. In order to usejt homoge¬ 
neous scale of v tur b, we adopted the calibration bv lKirbv et alj 
d2009b that provides v tm -b as a function of log g. No trend be¬ 
tween iwb and the K abundances are detected when the rela¬ 
tion bv lKirbv et al. ( 2009 ) is used. The K abundances are cor¬ 
rected for the departures from Local Thermodynamic Equi¬ 
librium (LTE) by applying non-LTE abundance corrections 
computed with the non-LTE radiative transfer code MULTI 
(modifie d version 2.3 ICarlssonl 19861) . MARCS model atmo¬ 
spheres dGustafsson et al.l 12008 ) and a new model atom of 


neutral potassium (Merle et al. in prep.). NLTE EWs are 
larger than in LTE, leading to NLTE K abundances lower by 
about -0.3 dex. 

Uncertainties in the derived abundances were computed by 
taking into account two different sources of errors: (i) the un¬ 
certainty arising from the EW measurement, as estimated by 
DAOSPEC, typically of the order of ~0.02-0.03 dex for GI¬ 
RAFFE targets and ~0.04-0.05 dex for UVES targets, and 
(ii) the uncertainties arising from the atmospheric parameters. 
For the latter source of error, we take into account the cor¬ 
relation among the parameters: we adopted 40 K as typical 
uncertainty in 7)n, as quoted by C06. A 7/n variation of ±40 
K changes the K abundance by ±0.06-0.07 dex and leads to 
a variation of log g smaller than ±0.1 (with a null impact on 
the K abundance). This log g variation leads to a change in 
Vturb of ±0.04 km/s, corresponding to a K variation of ±0.02- 
0.03 dex. Hence, we estimated a total uncertainty of about 
0.09-0.11 dex. 

4. RESULTS 

The mean difference between our radial velocity estimates 
and those by C06 is (AV r ) = -0.1 km/s with a standard de¬ 
viation of cta = 0.8 km/s. Three of the 119 targets have 
| AV r | >3.0 km/s, i.e. much larger than 3cta, hence they 
were flagged as possible binary stars and excluded from the 
following discussion (but shown in the plots as empty sym¬ 
bols). It is interesting to note that some of the candidate un¬ 
resolved binary systems identified here stand out as extreme 
outliers in the chemical abundance plots shown in Fig.[Q This 
suggests that some stars observed to lie out of the general 
trends between abundance ratio in GCs or other stellar sys¬ 
tems may not have a genuine anomalous composition but, in¬ 
stead, may have their abundances altered by their unrecog¬ 
nised binary nature. The abundance distribution of the re¬ 
maining 116 stars spans a sizeable range, from [K/Fe]=-0.16 
up to [K/Fe]=+0.33. We used the Maximum Likelihood al¬ 
gorithm described in Mu 12, that provides mean and intrinsic 
spread of a given abundance ratio by taking into account the 
uncertainties of each individual star. We found that the sample 
displays a statistically significant intrinsic star-to-star scatter 
in [K/Fe], cr^ /Fe] =0.05±0.01. 

Additional support for this conclusion comes from the be¬ 
havior of [K/Fe] as a function of [O/Fe], [Na/Fe], [Mg/Fe] and 
[Al/Fe] shown in Fig.[j]and Fig. [2] The GIRAFFE targets (for 
which only O and Na abundances are available) are shown as 
red circles and the UVES targets (for which O, Na, Mg and Al 
were measured) as green squares. Recently, ICarrettal (120141) 
re-analysed the stars already discussed in C09 and derived 
new Mg and Al abundances for the stars originally discussed 
in ICarretta et al] ( 2004 ) but not planned in our observations. 
However, two of these new stars are in common with our GI¬ 
RAFFE sample and we included their Mg and Al abundances 
(marked as filled triangles in Fig. 0, even if formally their 
analysis is based on a slightly different T t g- scale (with 7/n 
differences of 71 K and -18 K with respect to C06). 

All the four abundance ratios show clear and statistically 
significant correlations with the K abundance. In particu¬ 
lar, [Na/Fe] and [Al/Fe] correlate with [K/Fe], and [O/Fe] 
and [Mg/Fe] anti-correlate with [K/Fe]. Spearman corre¬ 
lation test gives probabilities that the variables are non- 
correlated of P(K-O) = 1.5 x KT 4 , P(K-Nci) = 3.9 x KT 4 , 
P(K-Mg) = 0.0665, and P(K-AI) = 0.0024. To better as¬ 
sess the significance of the correlations of K with Mg and 
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Al, we test the hypotheses that stars with [Mg/Fe]([Al/Fe])< 
0.0C+0.5) and those with [Mg/Fe]([Al/Fe])> 0.0(+0.5) (a) are 
extracted from the same parent distribution in [K/Fe], with 
a Kolmogorov-Smirnov test, and ( b ) are extracted from dis¬ 
tribution of [K/Fe] having the same mean, with Student’s 
and Welch’s tests. The hypothesis a) is rejected with P = 
1.8 x 10 -3 and hypothesis b) with P < 5.0 x 10 -5 . 

At the core of the Mg-K anti-correlation is the fact that 
the four Mg-deficient stars ([Mg/Fe]< 0.0) identified in 
NGC 2808, are also found to be K-enhanced with respect to 
Mg-normal stars. While the amplitude of the enhancement 
is relatively weak (~ +0.3 dex) the trend is analogous to that 
found in NGC 2419 by Mul2. Following the classification 
scheme set up by ICarrettal (120141) . Mg-deficient and K-rich 
stars should correspond to the Extreme population and Mg- 
normal and K-poor stars to the Primordial population. 

5. SANITY CHECKS 

We performed several sanity checks to assess the robustness 
to systematics of the results described above. 

Correlations with atmospheric parameters — Atmospheric 
parameters of stars along the RGB of a globular cluster are 
correlated. An increase of T e ff, coupled with the correspond¬ 
ing variations in log g and v tur b. produces an increase in all the 
abundance ratios. The effect of these changes on the abun¬ 
dance ratios is shown by the arrows in Fig. [T| and [2] The di¬ 
rections of the plotted vectors show that while a systematic in 
one of the atmospheric parameter can (in principle) play a role 
in the K-Na correlation (and, perhaps, also in the K-Al one), 
they are almost perpendicular to the directions of the K-Mg 
and K-O anti-correlations (see also the discussion in Mul2). 

Indeed, we noted that a weak residual trend between [K/Fe] 
and T e ff is present in our data (see panel (a) of Fig. [3}. How¬ 
ever, this trend with 7/if is clearly not at the origin of the 
K-O and K-Mg correlation, as it can be appreciated from 
Fig. 0 where we plotted with different symbols stars hotter 
and cooler than T e g= 4500 K. Both groups display the K-O 
correlation, with a small offset in K abundances fully con¬ 
sistent with the vector shown in Fig. Q] The application of a 
rigid offset of +0.08 dex in [K/Fe] to the stars of one group 
would lead to a full overlap of the two parallel sequences in 
the [K/Fe] vs. [O/Fe] plane. Note that such an offset corre¬ 
sponds to a variation of ^50 K, compatible with la uncer¬ 
tainty in T eff , but with a negligible impact on the abundances 
of other elements. A similar behavior occurs in the [K/Fe] vs. 
[Mg/Fe] plane, though the small sample prevents an interpre¬ 
tation as clean as that emerging from Fig. [3] The same effect 
is observed in the [K/Fe] vs. [Na/Fe] plane, explaining (at 
least partially) the large K abundance scatter measured at dif¬ 
ferent values of [Na/Fe]. It is especially relevant to note that 
the removal of this small systematic would made the observed 
correlations of K with other light elements even tighter than 
that shown in Fig. Q]and[2] Note that the three K-rich stars 
observed with UVES have a [K/Fe] scatter smaller than that 
measured among the GIRAFFE targets with similar [O/Fe] or 
[Na/Fe], This can be explained with the fact that these stars 
have similar 7j.fr, covering a range of about 100 K, while the 
entire sample of GIRAFFE targets covers about 800 K. 

NLTE corrections — [K/Fe] abundance ratios^ were re¬ 
derived using the NLTE corrections bv lTakeda et al.1 (I2002IR 
leading to a decrease of [K/Fe] of ~0.2 dex. The use of these 
NLTE corrections does not erase the observed correlations be- 

7 http://optik2.mtk.nao.ac.jp/fakeda/potassium_nonlte/ 


tween [K/Fe] and the other light elements, and the difference 
in K abundance between Mg-normal and Mg-deficient stars 
remains the same. The correlations are preserved also when a 
constant NLTE correction of -0.3 dex is adopted, as done by 
Mul2 in the case of NGC 2419. 

An alternative set of spectroscopic parameters — We de¬ 
rived spectroscopically the atmospheric parameters of the 
stars of NGC 2808 observed with UVES, thus repeating the 
analysis with a set of parameters fully independent of those 
by C09 that we used before. We combined the UVES spectra 
taken within our program, with those by C09 from the Red 
Arm 580 set-up, thus providing a total spectral coverage of 
about 4000 A . The atmospheric parameters were derived im¬ 
posing (a) no trend between Fe I abundances and excitation 
potential (to constrain T e ff), (b) no trend between Fe I abun¬ 
dances and EWs (to constrain v t urbX ( c ) the same abundance, 
within the uncertainties, from Fe I and Fe II lines (to con¬ 
strain log g). This new analysis, based on ~250 Fe I lines 
and ~ 20-25 Fe II lines, provides atmospheric parameters that 
match well those used above, with differences smaller than 
50-60 K in T e ff, ~ 0.1 in log g, and 0.2-0.3 km/s in v,. Conse¬ 
quently, the newly derived abundances are in good agreement 
with those presented above and the relevant trends with [K/Fe] 
remain essentially untouched. 

6. DISCUSSION 

The present analysis provides two main results: (1) the K 
abundance is not uniform among the stars of NGC 2808 but 
it exhibits a small but significant intrinsic spread; (2) K abun¬ 
dances are correlated with the abundances of the light ele¬ 
ments involved in the chemical anomalies ubiquitously de¬ 
tected in globular clusters (namely O, Na, Mg and Al). In 
particular, the four Mg-deficient cluster stars detected so far 
are all K-enhanced, implying the existence of a tight Mg-K 
anti-correlation similar to that observed in NGC 2419, albeit 
with a much smaller amplitude. This is the second case of a 
GC where a Mg-K anti-correlation is detected. 

As for other light elements in NGC 2808 (Carrettail2014l) 
the distribution of K is not unimodal, thus suggesting that 
the enrichment occurred in discrete events. The bi-modality 
of the K distribution is evident in the UVES sample. The 
GIRAFFE sample satisfie s all the criteria established by 
[Muratov & Gnedinl (120101) for a non-unimodal distribution 
(the so-called Gaussian mixture modeling, GMM, test). In 
particular, the GMM test rejects the hypothesis of unimodal 
distribution at the 99.5% confidence level. 

Stars with sub-solar [Mg/Fe] abundance ratios are quite 
rare in Galactic GCsj_ among the 19 GCs investigated 
by ICarretta et al] (120091) only NGC 2808 is found to 
harbor four stars with [Mg/Fe] <0; two were detected 
in M 54 (ICa rretta et aLl 120101), and four in Omega Cen- 
tauri dNorris & Da Costal 119951) . All these clusters (plus 
NGC 2419) are among the most massive of the entire 
GC system of the Milky Way and present the most ex¬ 
treme manifestations of the light element self-enrichment 
process typical of globulars. They display extended (anti- 
correlations between abundances of light elements, com¬ 
plex Horizontal Branch morphologies with extended blue 
tails, and are observed (or suspected) to host popula- 
tions with extreme He abundances (Y > 0.3, see, e.g. 


D’Antona et al. 2005; Dalessandro et al. 2008; Gratton et al. 

2011 

; Kina et al.l 2012; Dalessandro et al. 2011; Beccari et al. 

2013 

NGC 

, and references therein). Our results for NGC 2419 and 
2808 strongly suggest that, at least in these massive 
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Fig. 1.— [K/Fe] as a function of [O/Fe] (left panel) and [Na/Fe] (right panel). Red circles are the GIRAFFE targets and green squares the UVES targets. 
Empty circles are stars excluded from the analysis for their variable V r (likely binary stars). The arrows show the effects of a change by 40 K in T e ff and the 
corresponding variations in log g and v tU rb- 
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Fig. 2.— [K/Fe] as a function of [Mg/F e] (left panel) and [Al/Fe] (right panel). Green squares are the UVES targets, while the triangles are the additionally 
two stars in common with ICarrettal (2014 ). 


and metal-poor clusters, Mg-deficiency is associated with K- 
enhancement. Other clusters, where self-enrichment did not 
reach the stage leading to the production of Mg-deficient stars, 
probably were not able to synthesise enough Potassium to dis¬ 
play a detectable spread in [K/Fe] (ICarretta et al.ll2013l) . 

The evidence that the [K/Fe] correlates also with O, Na, 
and A1 supports the hypothesis that the spread in K is as- 
cribable to the same self-enrichment process that produces 
the observed abundance variations (and correlations) for these 
elements. Still, placing NGC 2419 in this scenario is not 
straightforward. In fact, while its K-Mg anti-correlation is 
consistent with being an extended version of that observed in 
NGC 2808, this is not the case for the K-Na and K-Al correla¬ 
tions, since, according to lCohen & Kirbvl(l2012h . at any given 
[Na/Fe] and [Al/Fe] there is a large spread of [K/Fe] (see also 
ICarretta et al.|[2013l for a thorough discussion). Perhaps this 


is somehow related to the different metallicity regime of the 
two clusters, whose [K/H] ratio for first-generation stars differ 
by ~ 1.0 dex (see also V 12 ). 

Within the scenario where the self-enrichment is driven 
by AGB stars, K production is expected to occur in the 
same hot bottom burning environment where the other chem¬ 
ical anomalies are produced (V12). In these models, the 
stars of the extreme populations of both NGC 2419 and 
NGC 2808 are thought to have formed directly from the pure 
ejecta of AGB/su per-AGB stars, without dilution with pris¬ 
tine g as (see e.g iD’Antona et all 120051: iDi Criscienzo et afl 
120111) . Potassium can be produced by proton cap¬ 
ture on Argon nuclei, through the thermonuclear chain 
36 Ar( p, 'y) il K(e + , v) 31 Cl(p, 7 ) 38 Ar( p, 7 ) 39 K, occurring at tem¬ 
peratures of about 10 8 K. Even if the chain able to produce K 
is known, several uncertainties affect our knowledge of this 
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Fig. 3.— Panel (a): trend of [K/Fe] abundance ratio with T e ff. The red 
line is a linear fit to the data. Panel (b): [K/Fe] as a function of [O/Fe]. Stars 
hotter and cooler than 4500 K are plotted with different symbols (red triangles 
and green circles, respectively) to highlight the effect of the trend between K 
abundance and 7^f on the [K/Fe]-[O/Fe] correlation. 

reaction. In particular, V12 had to increase the cross sec¬ 
tion of the reaction 38 Ar(p,y) 39 K by a factor of 100 to re¬ 
produce the K abundances measured in the Mg-deficient stars 
of NGC 2419. We refer the reader to V12 for a complete dis¬ 
cussion about the uncertainties in these reactions. 
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Note that no explicit predictions about K abundance anoma¬ 
lies are available for other models of self-enrichment (based 
on FRMS or binary stars) and we cannot use the observed 
scatter in the K content of NGC 2808 to discriminate among 
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depletion we observe in NGC 2808 stars (~ —1.3 dex). 
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